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ABSTRACT 

A spherically symmetric hydrodynamic stellar core collapse under gravity is time- 
dependent and may become unstable once disturbed. Subsequent nonlinear evolu- 
tions of such growth of hydrodynamic instabilities may lead to various physical con- 
sequences. Specifically for a homologously collapse of stellar core characterized by a 
polytropic exponent T — 4/3, we examine oscillations and/or instabilities of three- 
dimensional (3D) general polytropic perturbations. Being incompressible, the radial 
component of vorticity perturbation always grows unstably during the same homol- 
ogous core collapse. For compressible 3D perturbations, the polytropic index 7 of 
perturbations can differ from F = 4/3 of the general polytropic hydrodynamic back- 
ground flow, where the background specific entropy is conserved along streamlines and 
can vary in radius and time. Our model formulation here is more general than pre- 
vious ones. The Brunt— Vaisala buoyancy frequency J\f does not vanish, allowing the 
existence of internal gravity g~— modes and/or g+— modes, related to the sign of 7V^ 
respectively. Eigenvalues and eigenfunctions of various oscillatory and unstable per- 
turbation modes are computed given asymptotic boundary conditions. As studied in 
several specialized cases of Goldreich & Weber and of Cao & Fou, we further confirm 
that acoustic p— modes and surface f— modes remain stable in the current more gen- 
eral situations. In comparison, g~— modes and sufficiently high radial order g"*"— modes 
are unstable, leading to inevitable convective motions within the collapsing stellar in- 
terior; meanwhile sufficiently low radial order g"*"— modes remain stably trapped in 
the collapsing core. Unstable growths of 3D g— mode disturbances are governed dom- 
inantly by the angular momentum conservation and modified by the gas pressure 
restoring force. We note in particular that unstable temporal growths of 3D vorti- 
cal perturbations exist even when the specific entropy distribution becomes uniform 
and 7 = F = 4/3. Conceptually, unstable g— modes might bear conceivable physical 
consequences on supernova explosions, the initial kicks of nascent proto-neutron stars 
(PNSs) of as high as - 200 - 500 km s"^ up to > 1000 km s"^ and break-ups of 
the collapsing core, while unstable growths of vortical perturbations can lead to fast 
spins of compact objects, 3D vortical convections inside the collapsing core for possi- 
ble magnetohydrodynamic (MHD) dynamo actions on seed magnetic fields, and the 
generation of Rossby waves further stimulated by gravitational wave emissions. 

Key words: hydrodynamics — instabilities — waves — stars: neutron — stars: 
oscillations — supernovae: general 



1 INTRODUCTION 

Over past several decades, considerable theoretical and nu- 
merical efforts have been devoted to examine (in)stability 
properties of core collapsing stars before or around the emer- 
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gence of supernova (SN) rebound shock expanding inside 
a massive progenitor star (e.g. Goldreich & Weber 1980 - 
GW hereafter; Goldreich, Lai & Sahrling 1996; Lai 2000; Lai 
& Goldreich 2000; Blondin, Mezzacappa & DeMarino 2003; 
Murphy, Burrows & Heger 2004; Lou & Wang 2006, 2007; 
Burrows et al. 2006, 2007; Lou & Cao 2008; Cao & Lou 
2009, 2010; Lou & Wang 2011). A homologously collapsing 
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relativistically hot or degenerate gas sphere was first studied 
by GW to describe the pre-SN stellar core collapse process, 
in which the gas pressure p and the mass density p obey 
the conventional polytropic equation of state (EoS) p — Kp^ 
where F = 4/3 and k are two global constants. In physically 
realistic situations, both parameters F and k do not remain 
constant in space and time due to multiple complicated nu- 
clear processes including various neutrino productions, elec- 
tron captures etc. for the inner ~ I.SA/q iron core collapse 
according to numerical simulation results (e.g. Van Riper 
& Lattimer 1981; Burrows & Lattimer 1983; Bethe 1990; 
Hix et al. 2003). A decrease of the homologous stellar core 
mass may occur (e.g. Bethe 1990). One may then introduce 
a mean value < F > as a pressure weighted spatial average. 
Meanwhile, a few numerical simulations (e.g. Van Riper & 
Lattimer 1981; Burrows & Lattimer 1983, 1986; Bethe 1990) 
show that the specific entropy closely related to k does not 
vary much in time before the pre-shock stage for a pre-SN, 
and that the mean value < F > does not deviate from 
4/3 significantly. General relativistic simulations of Dim- 
melmeier et al. (2008) appear to confirm the near constancy 
of F before the emergence of a bounce shock. Moreover, other 
numerical simulations (e.g. Bruenn 1989a,b) show that the 
homologous solution for a stellar core collapse represents a 
fair first-order approximation. Based upon their homologous 
core collapse solution, GW further carried out a vorticity- 
free (i.e. a potential flow) three-dimensional (3D) perturba- 
tion analysis and claimed that no instabilities occur for such 
a dynamic stellar core collapse. Polytropic homologous core 
collapse solutions for F < 4/3 have been later explored by 
Yahil (1983), and 3D perturbation analysis on these gener- 
alized solutions later show that for F > 1.09 no instabili- 
ties arise (e.g. Lai 2000). We note that all these previous 
3D perturbation analyses are actually restricted to acous- 
tic p— modes as the Brunt— Vaisala buoyancy frequency J\f 
vanishes in these model considerations. 

Lou & Bai (2011) presented 3D isothermal perturba- 
tions in a self-similar isothermal nonlinear dynamic flow. 
Such 3D perturbations can be cast into a self-similar form 
to describe various possible 3D dynamic fiow structures and 
configurations. In fact, this isothermal requirement can be 
relaxed to allow a more general EoS. Lou & Cao (2008) 
generalized the self-similar or homologous core collapse so- 
lutions of F = 4/3 for a more general EoS by allowing a 
variable ft that is actually conserved along streamlines such 
that K becomes an arbitrary function of the independent 
self-similar variable which is a proper combination of radius 
r and time t; in other words, GW solutions only represent a 
special subclass for a conventional polytropic relativistically 
hot or degenerate gas of F = 4/3. Numerical simulations for 
the internal structure of massive progenitor stars and SN 
explosions (e.g. Bethe et al. 1979; Bruenn 1985; 1989a, b; 
Woosley, Langer & Weaver 1993; Woosley, Heger & Weaver 
2002) seem to indicate that the specific entropy distribution 
is generally non-uniform and can vary in space and time. 

The stellar iron core collapse, subsequent core bounce 
and shock emergence are complicated processes and involve 
several essential aspects of physics including various nuclear 
reactions, neutrino trapping and escape, energy transport, a 
sophisticated EoS (both at mass densities below and above 
normal nuclear density of ~ 2 — 3 x lO^^g cm~^), general rel- 
ativistic effects, variations of characteristic physical param- 



eters and so forth (e.g. Bethe 1990 and references therein). 
Hydrodynamic instabilities due to several relevant physical 
processes may appear before or during an SN explosion due 
to the core collapse of a massive progenitor star (e.g., with 
a mass M > IOM0). The development of numerical sim- 
ulation models and codes (mostly spherically symmetric) 
over five decades has reached various levels of sophistica- 
tion. Within the past decade in particular, the majority 
of simulation work has been multi-dimensional in nature. 
Recent multi-dimensional numerical simulations show that 
the growth of 2D and 3D disturbances can be significant 
and may be essential for the neutrino convection/heating 
and thus the viability of the subsequent SN explosion (e.g. 
Dessart et al. 2006; Buras et al. 2006a, b; Iwakami et al. 
2008; Nordhaus et al. 2010b). Among others, hydrodynamics 
(the so-called advective-acoustic "standing accretion shock 
instability" or SASfl see e.g. Foglizzo 2001; Blondin et al. 
2003; Blondin & Mezzacappa 2006) and neutrino-driven con- 
vections can lead to 2D instabilities. For example, it has been 
shown that the I — 1 mode is dominant and that a bipolar 
sloshing of the standing shock wave occurs, with strong pul- 
sational expansions and contractions along the symmetry 
axis. These numerical simulations indicate that the growth 
of initial perturbations seems to occur only after the nascent 
PNS is formed (e.g. Dessart et al. 2006; Nordhaus et al. 
2010b). In this context, contributions of the pre-bounce core 
collapse process to the growth of unstable perturbations 
should not be ignored, which could perhaps affect an SN ex- 
plosion directly or indirectly (e.g., by adjusting the process 
of neutrino convection and heating). In numerical simulation 
studies, it is not a trivial task to introduce self-consistent dy- 
namic fiuctuations during the process of core collapse phase; 
this is because in many situations, such perturbation solu- 
tions are not known or available a priori. Often, one has 
to rely on the inevitable imprecision of numerical codes as 
initial perturbations. This might miss unstable growth ef- 
fects of some perturbations during the earlier stellar core 
collapse phase. For example, Buras et al (2006a, b) studied 
partly effects of the pre-bounce perturbations and rotations 
by numerical simulations. They add arbitrary core collapse 
density perturbations, which would be more of acoustic na- 
ture, and find that the PNS convection starts slightly earlier, 
though nothing changes significantly. In the numerical sim- 
ulation of Iwakami et al. (2008) where small initial radial 
velocity perturbations are added, a linear growth phase be- 
fore the SN explosion is in fact observed, during which the 
amplitudes of the initial perturbations grow nearly expo- 
nentially with time t. Though they start their simulations 
from a post-bounce time, the stellar core is still contract- 
ing before the final explosion; this observation may still be 
traced back to implicate possible g— mode and/or vortical 
instabilities during the early stellar core collapse phase even 
when the shock has not yet emerged. It deems important 
to treat unstable perturbations with care during the early 
core collapse process. This paper is devoted to analyze the 



^ Physically, this SASI mechanism involves an amplifying feed- 
back cycle of entropy and vorticity perturbations, which are ad- 
vected inward and generate acoustic waves propagating outward 
to distort the shock front. Such acoustic waves grow between the 
PNS and accretion shock, or post-shock region. 
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growth of 3D unstable perturbations during the process of 
the stellar core collapse before the core bounce. 

In order to catch dynamic perturbation properties, we 
shall adopt in this paper a very much simplified yet still 
fairly challenging general polytropic EoS to describe a non- 
linear core collapse of spherical symmetry and perform 3D 
general polytropic perturbation (in)stability analysis. Bethe 
(1990 and extensive references therein) reviewed the simi- 
larity core collapse solutions (GW; Yaliil & Lattimer 1982; 
Yahil 1983) which appear to capture several gross features 
of numerical simulations. Though limited, such analytic and 
semi-analytic approaches may reveal certain physical as- 
pects more transparently. Likewise, we hope to derive use- 
ful hydrodynamic instability properties for a more general 
class of homologous core collapses and offer valuable phys- 
ical insights for further multi-dimensional numerical sim- 
ulation explorations, in spite of several simplifications and 
approximations involved in our hydrodynamic model formu- 
lation. Among others, such non-spherical effects have also 
been shown highly relevant to gravitational wave emissions 
during a violent dynamic stellar collapse involving rotating 
iron cores (e.g. Saen & Shapiro 1978; Dimmelmeier et al. 
2008 and extensive references therein). In numerical simula- 
tions, it is also valuable to introduce self-consistent pertur- 
bation solutions to initiate well-controlled code testing and 
to clearly identify physical effects. 

The 3D polytropic perturbation analysis for 7 = 4/3 on 
these generalized dynamic collapse solutions by Cao & Lou 
(2009) with a self-similar evolution/distribution of specific 
entropy show that qualitatively analogous to stellar oscil- 
lations about a static star (e.g. Cowling 1941; Cox 1976; 
Unno et al. 1979), all kinds of perturbation modes (i.e. 
acoustic p— modes, surface f— modes and internal gravity 
g— modes) exist, and instabilities can always develop for cer- 
tain g— modes of sufficiently high radial orders. Cao & Lou 
(2009) demonstrated that the sign of the Brunt— Vaisala 
buoyancy frequency squared JV^ , which determines the exis- 
tence of internal g— modes in static stellar models, still holds 
the key for g— modes during the self-similar dynamic evolu- 
tion phase of a collapsing stellar core. Prior perturbation 
analyses (e.g. GW; Lai 2000; Lai & Goldreich 2000) were 
based on the collapsing stellar core model with a globally 
constant specific entropy distribution, for which the gas flow 
is conventional polytropic and thus Af^ — 0, such that all 
g— modes are actually suppressed. Therefore, GW analysis 
only indicates that acoustic p— modes and f— modes remain 
stable in a homologous core collapse. Cao & Lou (2009) fur- 
ther confirmed the stability of p— modes and f— modes with 
a variable specific entropy evolution/distribution, but they 
clearly found that all g^ —modes and sufficiently high-order 
g"*"— modes are unstable during the core collapse of a massive 
progenitor star. They estimated and proposed that the in- 
stability of the lowest radial order I = 1 g— modes can lead 
to initial kicks of nascent PNSs with speeds of as high as 
~ 200-500 km s"^ up to > 1000 km s"^ (e.g., Scheck et al. 
2004, 2006, Nordhaus et al. 2010a, and references therein). 
Our PNS kick proposal differs from several other suggested 
mechanisms. For example, SN explosions in binary systems 
do not produce high enough kick speeds of PNSs (e.g. Lai 
et al. 2001). Hydrodynamic simulations for anisotropic mass 
ejection in an SN explosion have produced fairly low recoil 
velocities of NS (e.g. Janka & Miiller 1994) or started from 



the assumption that a dipolar asymmetry was present in 
the pre-coUapse iron core leading to a large anisotropy of 
an SN explosion (e.g. Burrows & Hayes 1996). Nevertheless, 
the physical origin of such pre-collpase fluctuations is not at 
all clear for the moment (e.g. Murphy et al. 2004). In fact, 
various kinds of hydrodynamic instabilities might be respon- 
sible for a large-scale deformation of the ejecta and globally 
aspherical SN explosions. Our unstable 1 = 1 g— modes dur- 
ing core collapse appear to be plausible for generating such 
dipolar asymmetry in envelope mass ejection. Note also that 
perturbation analysis of Chandrasekhar (1961) found the 
highest growth rates for the Z = 1, m = mode for thermal 
convection within a gravitating fluid sphere. 

Adopting the stellar collapse model with a constant 
specific entropy distribution as in GW, Cao & Lou (2010) 
considered non-isentropic perturbations with a perturbation 
polytropic index 7 / P = 4/3, for which J\f^ does not vanish. 
In a complementary perspective, the results of Cao & Lou 
(2010) also reveal instabilities of g~— modes and sufficiently 
high radial order g+— modes. 

To reveal behaviours of g— modes for the most general 
cases, we perform here a 3D perturbation analysis of the 
collapsing stellar core model by allowing not only the self- 
similar profile of the specific entropy to be arbitrary, but also 
the perturbation polytropic exponent 7 being different from 
P = 4/3 of the dynamic background. Under such more gen- 
eral conditions, the Brunt— Vaisala buoyancy frequency J\f 
does not vanish in general and g— modes exist as expected. 
The instability properties of g— modes are further explored 
in this paper, and several possible processes which may hap- 
pen during a pre-SN stellar core collapse are discussed in ref- 
erence to these instabilities, including the formation of PNSs 
and initial kicks of nascent PNSs. In addition, we also search 
for the mechanism of unstable g— modes in the core collaps- 
ing process of a massive progenitor star, and it is shown 
that the amplification of g— mode perturbations which are 
mainly convective are mostly due to the conservation of an- 
gular momentum. 

More importantly, 3D vorticity perturbations in a dy- 
namically collapsing stellar core with a uniform specific 
entropy distribution and perturbation polytropic exponent 
7 = r = 4/3 are in fact unstable. In the more general case 
of a variable specific entropy evolution/distribution and a 
perturbation polytropic exponent 7 P = 4/3, the radial 
component of vorticity perturbation always grows in am un- 
stable manner. The growth of such vorticity perturbations 
may lead to fast spins of collapsing stellar cores, vortical con- 
vections or circulations, as well as nonlinear Rossby waves 
over a spinning PNS and so forth. 

After an introduction of background information and 

our research motivation in Section 1, we present a hydrody- 
namic description for the dynamics of spherical self-similar 
general polytropic stellar core collapse in Section 2. For- 
malism for 3D general polytropic perturbations in a spheri- 
cally symmetric dynamic stellar collapse is shown in Section 
3. Model analysis and results are described in Section 4. 
Conclusions are summarized in Section 5. Details of spe- 
cific mathematical analysis can be found in Appendices A 
through D for the convenience of reference. 
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2 HOMOLOGOUS STELLAR CORE 
COLLAPSE OF POLYTROPIC GAS 

Fully nonlinear hydrodynamic partial differential equations 
(PDEs) include conservations of momentum and mass, the 
Poisson equation relating the gravitational potential "1? and 
the mass density p, and a general polytropic EoS (i.e. specific 
entropy conservation along streamlines); these PDEs are 



- + (u.V)u=--i^-V$ 
| + V.(pu) = 0, 



V^-l" = A-kGp 



+ u ■ V In 
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(2) 
(3) 
(4) 



where u(r, t), p{r, t), p{r, t) and $(r, t) are the bulk 
flow velocity, gas pressure, mass density and gravitational 
potential of a self-gravitating flow, respectively, the gradi- 
ent operator V is in terms of the 3D position vector r, 
G = 6.67 X fO~^ cm^ g~^s~^ is the gravitational constant, 
and r = 4/3 is the polytropic exponent for a dominant 
Fermi gas consisting of relativistic particles (mainly elec- 
trons, positrons, and neutrinofl in a radiation field). For 
a progenitor stellar core collapse prior to an SN explosion, 
the stellar core temperature is of order of several MeV or 
higher. This is indeed very hot in terms of gas temperature, 
i.e. ~ 10^° — lO'^^ K. Nevertheless, from the perspective of 
central electron chemical potential pe ~ 300 MeV and the 
neutrino chemical potential pv ~ 200 MeV, the core temper- 
ature of a few MeV is fairly low and the relativistic Fermi- 
Dirac gas remains highly degenerate (e.g. Bethe 1990). 

It is known that nonlinear hydrodynamic PDEs HJ — (HI) 
remain invariant under the following time reversal operation. 



t -t. 



P 



■ P, 



(5) 



such that a time-dependent hydrodynamic solution for a 
stellar core collapse may be reversed to effectively describe 
a dynamic outflow or expansion. 

In reference to Cao & Lou (2009), we first derive in the 
following self-similar dynamic collapse solutions of spherical 
symmetry and introduce the time-dependent spatial scale 
factor a{t) as 



ait) = p.(tr''^^ 



1/2 



(6) 



where the time-dependent pdt) and the constant coefficient 
Kc represent respectively the values of mass density p and 
specific entropy parameter k at the centre (hence a sub- 
script c) of a massive spherical progenitor star. With this 
spatial length scale factor a{t) defined by equation ([6]), the 

^ Neutrino interactions at sufficiently low densities (say a mass 
density p < few X lO^g cm~^) are rare enough that the stellar 
matter may be regarded as essentially transparent to neutrinos. 
For mass densities in excess of ~ 10l2-10"g 

cm ^ , the neutrino- 

baryon scattering becomes strong enough and neutrinos are effec- 
tively trapped, moving relative to the stellar matter only by dif- 
fusion. Neutrino-matter collisional cross sections are proportional 
to the square of the neutrino energy. 



vector radius r is replaced by a dimensionless vector radius 
X — r/a(t). From now on, we use uo(r, t), po{r, t), po{r, t) 
and ^o{r, t) to denote the dimensional unperturbed radial 
bulk fiow velocity, mass density, gas pressure and gravita- 
tional potential for the self-similar dynamic solution of the 
general polytropic flow background. The spherically sym- 
metric background flow takes the form of 



uo(r, f) = dx , 

Po(r, t) = ' a-' fix) 



po{r, t) = kpI^^ = 



<l>o(r, t) 



4/ < 



(7rG)2 

3 \ 1/2 



3 V vrG 



(7) 
(8) 

(9) 

(10) 



where K{r, t) bears the form of Kcg{x) with a constant Kc, 
and f{x) is a newly introduced reduced function of only x 
directly related the mass density po{r, t) of the dynamic 
background flow and is to be determined by hydrodynamic 
equations [see nonlinear ODE (|15() below]. Substituting ex- 
pressions (ff)) — H10|) for the dynamic background flow into 
nonlinear hydrodynamic PDEs ((T)) — ((J) by imposing spher- 
ical symmetry and without further approximations, we de- 
rive a set of coupled nonlinear ordinary differential equations 
(ODEs). It can be readily conflrmed that PDEs © and (g)) 
are automatically satisfled, and g(x) is allowed to take on 
a fairly arbitrary form, which gives the freedom to specify 
the self-similar evolution/distribution of the speciflc entropy. 
For the g{x) = 1 case, we would have the special subcase 
of a conventional polytropic relativistically hot or degener- 
ate gas with a global constant speciflc entropy parameter 
studied previously by GW. With a more general 
g{x) proflle for the speciflc entropy distribution/evolution, 
nonlinear PDE lU of momentum conservation leads to 



.gV'' ... 1 

— ^ I a a = — 

tin } X 



1 d 4 



4 d^/) 
3 da; 



(11) 



where the left-hand side (LHS) depends on time t only while 
the right-hand side (RHS) depends only on the independent 
self-similar variable a; = |x| = r/a{t). We should then re- 
quire both sides equal to a constant separation parameter 
4A/3, yielding the two separate nonlinear ODEs below 



da; 



TVG 

= Aa; 



1/2 



2.. 4A 



3 d 

4/3 da; 



(12) 



(13) 



ODE (|12p can be solved analytically to obtain an explicit 
temporal expression for the radial length scale factor a{t) 
necessary for deflning the independent self-similar variable 
a: = |x| = r/a(t) by equation ([6} and the radial flow uo = dx 
of the dynamic background by equation ((TJ . The solution for 
the spatial scale factor a{t) from ODE (|12|l for a homologous 
stellar core collapse (GW) is simply 



a(t) = (6A) 



1/3 



1/6 



,2/3 



(14) 
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where a positive separation constant A > is a physical 
requirement. Substituting ODE (|13p into Poisson equation 
((3| that relates the gravitational potential $ and the mass 
density p, we readily attain a nonlinear ODE of f{x) below 



_l__d_ 

dx 



p Ax 



(gf 



Af = 4A . 



(15) 



We can prescribe substantially diflFerent yet plausible forms 
of g{x) for specific entropy evolution/distribution and nu- 
merically integrate nonlinear ODE (llSp with proper bound- 
ary conditions using the standard fourth-order Runge-Kutta 
integration scheme (e.g. Press et al. 1986). For instance, it is 
possible to impose f{x},) = at the moving boundary x = xt 
for a homologously collapsing stellar core (e.g. a PNS). Once 
f{x) is solved satisfying the specified boundary conditions, 
the mass density po{r, t) and other physical variables char- 
acterizing the general polytropic dynamic background flow 
can all be derived accordingly [see eqns ((7)) — (ITO)) and ([13}]. 

For the A = case, nonlinear ODE (|15p reduces to the 
general polytropic Lane-Emden equation (e.g. Eddington 
1926; Chandrasekhar 1939; Cao & Lou 2009). As an addi- 
tional necessary check of the consistency, ODE psp reduces 
to equation (16) of GW precisely as expected for g{x) — 1. 
We would set /(O) = 1 as a consistent normalization for 
nonlinear ODE psp . In general, Cao & Lou (2009) require 
g'{0) = and /'(O) — for a finite sound speed at the cen- 
tre, where the prime ' indicates the first derivative in terms 
of the independent self-similar variable x. As first shown by 
GW and confirmed by Lou & Cao (2008), there exists a 
maximum acceptable value of A > 0, denoted by Am, for 
f{x) to vanish at a finite value of a:: = xt; in the parameter 
regime of A > Am, f{x) would have a minimum but will not 
reach zero within the semi-infinite interval ^ x < +oo. By 
numerical computations for the special case of g{x) — 1, this 
range for sensible A values is ^ A Am = 0.00654376. 

In our scenario, we identify this Xt as the contracting 
surface of the inner homologously collapsing stellar core and 
the enclosed core mass within xj, remains unchanged. In ear- 
lier numerical simulations on the stellar iron core collapse 
(e.g. van Riper & Lattimer 1981; Burrows & Lattimer 1983; 
Dimmelmeier et al. 2008), both F and k change in a com- 
plicated manner along mass parcels as dictated by several 
nuclear processes. One consequence of variable k and F in 
a collapsing star is that the mass of the inner homologously 
collapsing core decreases during the stellar core collapse to 
form a PNS with a mass of no more than ~ 0.5Mq at the 
end of collapse. In our dynamic model of stellar core collapse, 
F = 4/3 does not change and K(r, t) — Kcg{x) evolves in a 
self-similar yet fairly arbitrary manner. As the enclosed mass 
M within Xb remains constant, we cannot directly model the 
gradual mass reduction of the inner homologously collaps- 
ing core. By adjusting parameters of our model, we could 
describe the evolution of inner homologously collapsing core 
with low, intermediate and high masses, separately. We may 
also relax the requirement of f{xi,) = at x = xt and choose 
instead f{xb) — /min > for A > Am, where /min > is the 
minimum value of f{x) > 0. As such, this leaves freedom 
to match the inner homologously collapsing core with the 
outer envelope collapse. 

Given the idealizations and limitations of our model 
formalism, we invoke such a self-similar background flow so- 
lution to represent a inner stellar core collapse (GW; Lou & 



Cao 2008; Cao & Lou 2009, 2010). Our motivation here is 
to examine oscillations and instabilities of 3D general poly- 
tropic perturbations during such a hydrodynamic core col- 
lapse of spherical symmetry. The results of our model anal- 
ysis would reveal interesting physical effects to be further 
explored by multi-dimensional numerical simulation analy- 
sis despite of the relative simplicity of our model scenario. 



3 GENERAL 3D DISTURBANCES IN A 
POLYTROPIC PROCESS 

We now introduce general polytropic 3D perturbations to 
the dynamic background flow of self-similar core collapse de- 
scribed and summarized in Section 2. Specifically, we adopt 
the following forms for general polytropic 3D perturbations, 



u = a(x + vi) , 
P = Po(l + /i) , 

3 ^ 1/2 



$ : 



3 V ttG 



(16) 
(17) 
(18) 

(19) 



where po{r, t) and poij', t) of the dynamic background flow 
are expressed by eqns ([8]) and ((9|, respectively. In the back- 
ground plus perturbation expressions (|16p — (|19p . four di- 
mensionless variables with subscript 1 viz. vi, /i, Pi, ip\ 
(together with the pertinent multiplicative factors in front) 
represent 3D perturbations superposed onto the correspond- 
ing hydrodynamic background variables uo, po, Po, ^o, re- 
spectively. For a more general consideration, we allow such 
3D perturbations to be a general polytropic process with a 
polytropic exponent 7 being different from the polytropic 
exponent F = 4/3 characterizing the spherically symmetric 
dynamic background for a stellar core collapse. If 7 actu- 
ally takes on the value for the ratio of specific heats, then 
such 3D perturbations would be adiabatic. If7 = F = 4/3 
and with a variable g{x), we come back to the perturbation 
analysis of Cao & Lou (2009). 

Hereafter, the modified nabla operator V is with re- 
spect to the spherical polar coordinates {x, 0, ip) instead 
of (r, 9, ip) unless otherwise stated. After substituting ex- 
pressions die])- (US]) into PDEs (Il|-(|4]) with the standard 
linearization procedure and defining a logarithmic time 



r = ~ln\t\ 



(20) 



to consistently remove apparent t dependence in the coeffi- 
cients of the self-similar hydrodynamic background in per- 
turbation equations, we readily obtain linearized equations 
in the following forms of 

dvi vi _ 1 

3" " 4A/3 



V(5/"/Ji) - fiS/igf) 



3A 



dh 
dr 

dr 



3f /i 



2V ■ (f vi) 



3/3 



7(V ■ vi) + vi ■ 4 



V/ 



g 



(21) 
(22) 
(23) 

(24) 
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for 3D general polytropic perturbations in vector momen- 
tum equation, Poisson equatior(f|, mass conservation, and 
the specific entropy conservation along streamlines, respec- 
tively. In deriving equation (|24[) . we have changed F in PDE 
Q into 7 for 3D perturbations, as mentioned in the pre- 
ceding paragraph. We have also arranged the sign of r in 
definition pO|) such that as t — >■ —0^, r approaches positive 
infinity +oo. Taking the partial derivative with respect to r 
and the curl V x operation of momentum perturbation equa- 
tion ((2T|) . and substituting linear PDEs ((23)) and ((24)) into 
this resulting equation, we immediately arrive at a vector 
equation in terms of vi without involving the other three 
perturbation variables j3\, fi and tpi, namely 



2A 3 



d_ 

dr 



(V X vi) 



(4 - 3-y)gVf + fVg 



V(V ■ vi) + 3Vg X V(vi ■ V/) - 3V/ x V(vi • Vg) . (25) 

For the separation of variables in PDE ([25}, we may write 
vi = g(r)vi where g(r) is a function of only t and vi 
depends on {x, 0, (f}}. It follows that 

4m r 

— (V X vi) = [(4 - 37)pV/ + fVg x 

V(V ■ vi) + 3V3 X V(vi ■ V/) - 3V/ x V(vi • Vg) (26) 
and 



3^ - 1^: I^W 



dr 



2m 



(27) 



where m is the separation constant. For g(r) oc exp(?7r) in 
ODE ((27)), we obtain Zrf ~ -q ~ 2m/(3A) = and thus two 
roots ■n=[l±{l + 8m/A)^/^]/6. The RHS of equation ((25)) 
does not vanish in general, except for the special subcase 
when 7 = F = 4/3 and g{x) = 1 (i.e. V(? = 0) happen simul- 
taneously. In that special subcase, the RHS of equation (|25|) 
vanishes and the vorticity perturbation (V x vi) has an un- 
stable temporal growth e^^'^ during the dynamic stellar core 
collapse, although acoustic p— modes and f— modes remain 
stable (see GW and Cao & Lou 2009). We show presently 
with more details in Section 4.1 that dimensionless vortic- 
ity perturbation (V x vi) unstably grows in time as 
with a fairly arbitrary dependence on (a;, Q, (j>). Physically, 
such growing vortical motions satisfying the mass conserva- 
tion (see Appendix A) may give rise to rapid core rotations 
and 3D circulations or convections inside a self-similar col- 
lapsing stellar core. This can provide a sensible mechanism 
of amplifying seed magnetic fields by magnetohydrodynamic 
(MHD) dynamo actions in a collapsing stellar core. 

We solve the 3D perturbation problem in the modified 
spherical polar coordinates (x, 6, if) and presume a suffi- 
ciently general form of flow velocity perturbation vi as 



Vl = BrWr + V± (xiUt) -|- V X (lUrotBr) 



(28) 



where er stands for the unit radial vector, Wr, wt and Wr ot 
are three dependent variables characterizing the flow veloc- 
ity perturbation vi; and the modified transverse gradient 
operator is defined by 

^ In reference to equation I I22II , it appears that perturbed Poisson 
equation (25) of GW contains a sign typo (see correct equation 
27 of Cao & Lou 2009). 



Vx=e,-^+e^^-|- , (29) 
xoO X sm dip 

where eg and represent two orthogonal unit vectors in 
the 6 and (p directions, respectively. In expression ((28)), it is 
easily seen that physically Wr stands for the radial velocity 
perturbation, wt governs the tangential (i.e transverse to 
the radial direction) potential flow perturbation, and that 
Wrot is responsible for the tangential vortical (non-potential 
part) flow perturbation. By the superposition principle, we 
decompose the angular dependence of Wrot into the spherical 
harmonic function Yim{9, ip), characterized by a pair of two 
integral indices / and m. It follows from expression (|28() that 

er ■ (V X Vl) = -Vlwrot = l{l + l)Wrot / X^ . (30) 

It is then clear that for nonzero radial vorticity perturba- 
tion, we must require 1^1, also expected on the intuitive 
ground. It can be readily seen that the RHS of the radial 
component of the curl of momentum perturbation equation 
(|21|l vanishes, yielding an independent differential equation 
of Wrot below in reference to expression (|30() . 



3^- 
or 



1 VxWrot = 



3^-1 ]l{l + l)Wrot/x' =0 .{31) 



for each spherical harmonic component Yiin{9, ffi)- Solving 
this radial vorticity perturbation equation of Wrot, a grow- 
ing time dependence e^'^^ = for all spherical har- 
monic components of Wrot can be readily derived. This tem- 
poral factor naturally shows an increase of the amplitude as 
t — > —0" for the radial component of vorticity perturbation 
by expression (|30p . A non- vanishing Wrot leads to a per- 
turbation solution for vortical motions with fairly arbitrary 
angular dependence in {6, (j>) superposed onto a dynamically 
collapsing stellar core without involving the radial compo- 
nent of velocity perturbation. These perturbation modes of 
Wrot physically correspond to various vortical motions trans- 
verse to the radial direction Gr. We emphasize that the only 
constraint on Wrot is ODE (|3ip . with no other spatial con- 
straints (in terms of x) on Wrot- Therefore perturbation flow 
component patterns of Wrot with various combinations of 
{I, m} are all allowed and can be superposed together with 
a set of proper coefflcients to form fairly general or arbi- 
trary circulation patterns. Meanwhile, it is quite clear that 
all {I, m} vortical modes share the same temporal growth 
factor e^''"^, and they may form rather random tangential 
vortical flow patterns according to the superposition princi- 
ple. For the particular case of / = 1, m = as an example, 
we take Wrot = x^Yifi{6, ip) = x^cosO, and then equation 
(|28|l gives an expression for the corresponding velocity per- 
turbation as Vl — xsinO e^; and this actually describes a 
global spin motion to perturb the radially collapsing stel- 
lar core. In the top panel of Figure 1, we show a tangential 
vortical flow pattern in which both a global spin of the stel- 
lar core and a relatively isolated vortex coexist. It is easy 
to understand the time dependence e'^^^ here by referring 
to expression (|16p . As can be seen, this type of vortical ve- 
locity perturbation is given by ui = avi oc e'^^^^ oc . 
Noting also that Wrot represents generally vortical pertur- 
bations (e.g. the above perturbation example of rotational 
motion) for which the angular momentum is conserved, it is 
then clear that the amplitude increase of this kind of vortical 
velocity perturbation results from the conservation of angu- 
lar momentum and the rapid radial contraction of the stellar 
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Figure 1. Streamlines of transverse circulations for three illus- 
trating examples of tangential vortical motions over a spherical 
surface. The patterns of such vortical perturbation motions over 
spherical surfaces can be fairly arbitrary, while the magnitude 
of velocity perturbation increases inversely proportional to the 
radius of the collapsing stellar core due to the conservation of 
angular momentum. In the top panel, both a global spin of the 
collapsing stellar core with an embedded vortex are displayed. 
For the case of Wrot oc Y2^i{0, ip), the circulation pattern over a 
spherical surface is shown in the middle panel; and that for the 
case of Wrot oc 15,1(6, (p) is displayed in the bottom panel. 



inner core size. For such growing vortical perturbation so- 
lutions, the perturbations in stellar pressure, mass density 
and gravitational potential (/i, /3i and ijji) all remain zero 
(i.e. this is a type of incompressible flow perturbations). 

The temporal growth factor e^-^'^ here is faster than the 
temporal amplification factor e^^^ for the oscillatory com- 
pressive normal modes [e.g. acoustic p— modes with 7 = T = 
4/3 and g{x) — 1]; this amplitude increase for all oscillatory 
p— modes corresponds to an adiabatic amplification of sound 
waves due to compression of the core collapse according to 
GW. In contrast, this faster e'^^^ growth factor reveals a 
non-oscillatory unstable amplitude increase of radial vortic- 
ity perturbation during the self-similar dynamic phase of a 
stellar core collapse. By our model scenario, the fast spin of 
a nascent PNS or a stellar mass black hole (identified as the 
collapsed stellar core remnant) is most likely triggered and 
gained during such an unstable growth of i = 1, m = vor- 
tical perturbations about the core collapse and concurrent 
nonlinear dynamic interactions. 

There are alternative mechanisms proposed by others 
which might also contribute to the fast spin of a nascent 
PNS. Several years ago, Blondin & Mezzacappa (2007) sug- 
gested that the 3D growth of Z = 1 spiral modes (see their 
figures 1 and 2) of the SASI during the post-bounce delay 
phase may significantly spin up the core matter behind the 
stalled shock. We note more specifically that this type of 
SASI fluctuations is mainly of acoustic nature. In their 3D 
simulations, a nascent PNS has gained a considerable spin 
velocity with final periods of ~ 70ms for moderately rotat- 
ing progenitors and ~ 50ms for non-rotating progenitors. In 
contrast, more recent numerical simulations (e.g. Iwakami 
et al. 2008; Rantsiou et al. 2011) seem to show that such 
I = 1 spiral mode mechanism does not work so effectively 
and may only induce possible PNS spin periods of seconds. 
Blondin & Mezzacappa (2007) invoked several simplifica- 
tions in their 3D simulations; they found that the nonlin- 
ear evolution of the SASI is dominated by low-order non- 
axisymmetric mode characterized by a spiral flow pattern 
beneath the accretion shock during the post-bounce delay 
phase (e.g. Blondin & Shaw 2007; Fernandez 2010). Strictly 
speaking, they did not follow collapse itself, but started in a 
steady-state post-bounce configuration. Also, neutrino heat- 
ing and cooling were neglected and the fiow was assumed to 
be isentropic in Blondin & Mezzacappa (2007). In short, the 
issue of fast PNS spin remains open and should be explored 
further for vortical instabilities proposed here. 

Conceptually for the dynamic core collapse of an ex- 
tremely massive progenitor star (e.g. Lou & Wang 2006, 
2007; Hu & Lou 2009; Lou & Wang 2011), a fast spinning 
central black hole may similarly form resulting from an un- 
stable growth of such I = 1 and m = vortical perturba- 
tions. As an example of estimation, for a core collapse in 
which the stellar radius contracts by a factor of ~ 10~^, 
an initial stellar spin angular velocity of ~ 1 — 10 cycles per 
day would lead to a spin angular velocity of ~ 10 — 100 cycle 
per second after the core collapse. Such effects might have 
already been observed in some numerical simulation works, 
e.g. Buras et al. (2006b), in which an initially specified stel- 
lar angular velocity is significantly amplified. This effect may 
considerably or even dominantly spin up a nascent PNS or 
a stellar mass black hole among other alternatives. More- 
over, circulation speeds of any small or large vortices tan- 



8 Yu-Qing Lou & Biao Lian 



gent to spherical surfaces within a dynamically collapsing 
stellar core will be amplified according to this law; in 
other words, cyclones on spherical layer surfaces would spin 
faster and faster during such a dynamic core collapse. 

It is conceivable that nonlinear interactions among such 
large-scale vortices within spherical layers and stellar core 
rotation would generate and sustain Rossby-type inertial 
waves (e.g. Haurwitz 1940; Lou 1987, 2000, 2001). Gravi- 
tational wave radiation carries angular momentum and may 
further excite stellar Rossby waves to set an upper limit on 
the spin rate of young neutron stars (e.g. Andersson et al. 
1999; Dimmelmeier et al. 2008). If the stellar core collapse 
process is accompanied with seed magnetic fields, these fast 
increasing vortex perturbations nonlinear interactions might 
also induce magnetohydrodynamic (MHD) dynamo actions 
to enhance magnetic field strengths (e.g. Thompson & Dun- 
can 1993; Hu & Lou 2009). 

In the following general polytropic 3D perturbation 
analysis, we thus exclude the part involving Wrot, as we 
have already shown that it leads to unstable vortical per- 
turbation solutions (i.e. radial component of vorticity per- 
turbation) including the growth of rotational instability in a 
homologously collapsing stellar core. To separate the tempo- 
ral component from the angular components of perturbation 
variables, we can consistently presume that j3\, fi, Wr 
and wt bear the two-dimensional (2D) angular dependence 
of spherical harmonics Yi^{8, tp) and the temporal depen- 
dence e"^ with n being an additional parameter to be de- 
termined. Combining linearized equations (|23p . 1)24^ and the 
two angular components of momentum perturbation equa- 
tion (|21|) to eliminate the two dependent variables /3i and /i , 
it is straightforward to derive the following linear ODEs for 
3D general polytropic compressible perturbations, namely 



1 d% , , 1(1 + 1) , , . „,2 



— — a; Wr) 

dx 



lil + l)f^ 



-Wt 



3/Vv 



, 



-mxwt - ^gf 



4 

-m 

3 



Wr —{xWt) 

dx 



1 d 2 ^ 1(1+ 1) 

— — (a; Wr) Wt 

x^ dx X 



i^gf + f9')wr~2{ni>i)=0 
+ [(4 - 37)<?,f ' + fg'] 



1 d 2 ^ ^(^+1) 

— — (a; Wr) 

x-' dx X 



Wt 







(32) 



(33) 



(34) 



In three linear perturbation ODEs (|32|) — (|34|) above, /' and 
g' are respectively the first derivatives of f{x) and g{x) with 
respect to x, and m — 3An(3n — l)/2 is the eigenvalue pa- 
rameter to be determined by proper physical requirements. 
This definition of the eigenvalue m agrees with that of Cao 
& Lou (2009). In fact, equation (|34|) is simply the angu- 
lar component of vorticity perturbation equation (|25l) . Note 
that we can actually regard the product mpi as one depen- 
dent perturbation variableQ Perturbation ODEs (|32l) — H34p 
are equivalent to vector equation (|26|) . 



For physically sensible perturbations, the relevant 
boundary conditions take the following form of 

nt/ii oc x\ Wr = Iwt oc s'^^ for x 0^ , 



for X 
for X 



nipi (X X ^ ' 
mxwt , n%lJi 

— gf Wr ^ for x = Xb; (35) 

here stands for the value of x at the moving boundary 
of the collapsing stellar core where the mass density po van- 
ishes. These boundary conditions are imposed to avoid sin- 
gularity in perturbation solutions as a:: — >■ O""" and x — ^ -|-oo. 
The condition at the contracting boundary Xb requires a zero 
Lagrangian pressure perturbation there. 



4 MODEL ANALYSIS AND RESULTS FOR 
THREE-DIMENSIONAL PERTURBATIONS 

4.1 A re-visit to GW model analysis and results 

We first re-visit the special conditions adopted by GW, in 
which 7 = r = 4/3 and g{x) = 1 for a constant specific 
entropy. It follows that our equation (|34|) reduces to 



= . 



(36) 



For this reason, GW choose to set Wr ~ d{xwt) /dx with the 
eigenvalue parameter m 7^ 0, leading to a kind of 3D po- 
tential flow velocity perturbation solutions representing only 
acoustic p— mode and f— mode oscillations in a homologously 
collapsing stellar core. This is because internal g— mode os- 
cillations necessarily involve vorticity perturbations. Appar- 
ently, vorticity perturbations were completely excluded in 
the model analysis of GW. 

In fact, there exist other perturbation solutions by set- 
ting m = in equation (|36[) . corresponding to either n — 
or n — 1/3 yet leaving the vorticity perturbation factor 
[wr — d{xwt)/dx] being fairly arbitrary in terms of x [see 
also equation (|25|) with a vanishing RHS]. In the former 
case of n = 0, all perturbations are neutrally stable, while 
the pressure, density and gravitational potential perturba- 
tions do not need to vanish in general. We now focus on the 
latter case of n = 1/3, because in this case the perturbation 
is amplified during the self-similar core collapse process and 
is faster than the amplitude increase e^'^® for the oscillatory 
acoustic p— modes due to adiabatic amplification of sound 
waves by compression in the core collapse. It can be demon- 
strated (see Appendix A for details) that to make equations 
H32|l and (I33|l meet this demand for n — 1/3, we must have 



— — (a: Wr) ^ —Wt + 3f Wr = 

x-' dx X 



(37) 



corresponding to mass conservation without mass density 
perturbation, i.e. the divergence of velocity perturbation 
vanishes. This can be easily confirmed by substituting equa- 
tion (|37|l into perturbation ODEs (|32p and (|33p . For these in- 
compressible vorticity perturbation solutions, the pressure, 
density and gravitational potential are not disturbed. How- 
ever, the time dependence of e^^^ for n — 1/3, which is faster 



* The dependent perturbation variables Wr and wt here are pro- 
portional to Wr and wt respectively of Cao & Lou (2009). For 



precise notation definitions, we should identify our Wr and wt 
here with — 3nuir/2 and — 3nuit/2 of Cao & Lou (2009). 
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Figure 2. Eigenvalues of m for the case of A = 0.001, 7 = 1.37 > 
r = 4/3 and g{x) given by expression II39I I as an example of illus- 
tration. The moving boundary of the stellar core is numerically 
determined at x;, = 6.9561. In this case, both g+— modes (only 
three branches are shown) and g~— modes (only one branch is 
shown) exist. The eigenvalues of g~— modes (to the right of the 
vertical dash-dotted line) are quite close to zero. The eigenvalues 
of the first two orders of p— modes are on the left side. There are 
more branches of p— modes of higher orders with smaller eigenval- 
ues of m not shown in this figure. The unique f— mode branch with 
I ^ 2 separates p— modes and g— modes (n.b. the I = 1 f— mode 
does not exist). If 7 value is higher, eigenvalues of g~— modes 
tend to vanish. 

than the contraction amphfication factor e^''^ for oscillatory 
acoustic p— modes in the dynamic background core collapse 
(see GW), indicates a relationship between these growing 
vorticity perturbation solutions and the angular momentum 
conservation (see the discussion on the temporal factor e^''^ 
in Section 3). In fact, these perturbation solutions represent 
circulatory motions with non-vanishing vorticities, so it is 
natural that they develop according to the angular momen- 
tum conservation. In summary, we are left with ODE (|37|) 
where [wr — d{xwt)/dx] is an arbitrary function of x. By 
equation (|25|) . this describes the growth of vortical convec- 
tive motions, including in particular components of vorticity 
perturbation transverse to the radial direction. 

We note that the mass density is not significantly dis- 
turbed in convective motions or vortex circulations; this fea- 
ture also holds approximately for internal g-modes to be 
discussed further. In Buras et al. (2006b), they introduced 
arbitrary initial density fiuctuations in the pre-bounce core 
collapse to examine effects of pre-bounce perturbations, and 
they found nothing obviously different except that the PNS 
convection starts slightly earlier. According to our analysis 
above, their perturbations are more likely of acoustic na- 
ture rather than convective motions and grow much slower, 
having relatively minor effects on the post-bounce phase. 

4.2 Computational Analysis and Explorations 

In order to solve the perturbation eigenvalue problem posed 
here, we implement the procedure outlined below: first we 
use the fourth-order Runge-Kutta method (e.g. Press et al. 
1986) to numerically integrate equation (|15p for obtaining 
the self-similar background function f{x) which is directly 
related to the dynamic background mass density po{r, t) by 



Figure 3. Eigenvalues of m for the case of A = 0.001, 7 = 
1.30 < r = 4/3 and g(x) given by expression II39I I. The contract- 
ing boundary of the stellar core is numerically determined to be 
at a:(, = 6.9561. Compared to those shown in Figure 2, the eigen- 
values of m in this case are shifted towards the positive direction. 
The vertical dash-dotted line at m = separates nearby branches 
of g"*"— modes (left; only three branches are shown) from nearby 
branches of g~— modes (right; only three branches are shown). 
When the perturbation polytropic index 7 becomes sufficiently 
small, eigenvalues of g"*"— modes would completely disappear. 

equation ((SJ. With this ready, we systematically discretize 
linear ODEs (|32p — (I34|l with a proper mesh to cast our 3D 
perturbation problem to a matrix eigenvalue problem with 
specified boundary conditions H35p . The inverse iteration 
method (e.g. Wilkinson 1965) is then employed to search for 
the eigenvalues of m and the corresponding eigenfunctions 
(see Appendix C for details of the numerical scheme). 

We first solve the cases of globally constant specific en- 
tropy [i.e. g{x) = 1] in space and time with 7 = F = 4/3, 
and the results of analysis agree quite well with those of 
Cao & Lou (2009) again confirming the errors of p— mode 
eigenvalues as computed by GW. We have also checked and 
confirmed the results of the sample solutions in Cao & Lou 
(2009) for several variable forms of g{x) and 7 = F = 4/3 
and in Cao & Lou (2010) for a constant g(x) = 1 and several 
7 7^ F = 4/3. With the assurance of a few necessary consis- 
tencies, we continue to apply our numerical codes to further 
explore the more general cases with a nontrivial variable 
g{x) 7^ 1 and the polytropic exponent 7 7^ F = 4/3. 

In stellar oscillations about static stars, there exist two 
types of g— modes (e.g. Cox 1976) which are distinguished 
by the sign of the square of the Brunt— Vaisala buoyancy 
frequency defined by 

/9^_19^\ 
\ or 7 dr / 

where po and po are the background gas mass density and 
thermal pressure, respectively, C/ > is the local gravita- 
tional acceleration and 7 is the polytropic exponent of 3D 
perturbations. For stellar oscillations, if M'^ is positive ev- 
erywhere, only g^— modes occur, whereas if M'^ is negative 
everywhere, only g~— modes can be found. 

For our self-similar core collapse case under consider- 
ation, definition (I38p for A/"^ can be extended to allow po 
and Po being also time-dependent for a dynamic background 
fiow. Our numerical exploration demonstrates that if M'^ is 
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Figure 4. The eigenfunctions for the lowest order g+— mode with 
A = 0.001, 7 = 1.37, 1 = 2 and g{x) given by expression II39I I (see 
Fig. 2 for sample eigenvalues) . The moving boundary of the col- 
lapsing stellar core is at Xf, = 6.9561. The oscillation is primarily 
confined within the central region of the collapsing stellar core. 



Figure 5. A higher order mode with A = 0.001, 7 = 1.37, I = 
2 and g{x) given by expression II39I I. The moving boundary of the 
collapsing stellar core is at xi, = 6.9561. The oscillation overcomes 
the barrier where g'{x) < and approaches the outer part of the 
stellar core. This appears to be a kind of wave tunneling effect. 



positive everywhere, the eigenvalues m are always negative 
and only g"''— modes occur, whereas if Af^ is negative every- 
where, the eigenvalues m become always positive and only 
g~— modes can be found. As an example of illustration, we 
choose the dimensionless function g{x) for the specific en- 
tropy evolution/distribution in the form of 

g{x) = 1 + O.lx"^ exp(-a;V2) , (39) 

and take A = 0.001 for the self-similar hydrodynamic back- 
ground. In Figures 2 and 3, we show the eigenvalues m for 
different spherical harmonic degree / computed by setting 
7 = 1.37 > 4/3 (Fig. 1) and 7 = 1.30 < 4/3 (Fig. 2), respec- 
tively. In our numerical exploration, it could be easily seen 
that the eigenvalues m are shifted towards the positive di- 
rection when 7 decreases gradually. In Figures 2 and 3, both 
g"*"— modes and g~— modes exist. However, when the value 
of 7 is high enough (low enough), eigenvalues of g~— modes 
(g'*'— modes) will disappear by examining the trends of ex- 
tensive numerical results. This can also be seen from defini- 
tion (|38|) of the Brunt— Vaisala buoyancy frequency squared 
J^^ . Note that the eigenvalues m = 3An(3n — l)/2 of 
g— modes are all fairly close to 0, corresponding to approx- 
imately temporal growth factors around e^^^ = for 
the amplitude of oscillations when n ~ 1/3. 

In Figure 4, we show the lowest order of g"'"— modes for 
the case 7 = 1.37. The solid curve stands for the scaled 
eigenfunction xf^^'^Wr (related to the radial velocity per- 
turbation), and the dash-dotted curve represents the scaled 
eigenfunction x^f^^'^Wt (related to the transverse velocity 
perturbation). We also show the curve of the function g'{x) 
by the dotted curve [see expression (|39|) for the form of g{x)] . 
It can be seen that the oscillatory wave pattern of the eigen- 
functions are mainly trapped within the collapsing stellar 
core. However, the patterns are no longer confined inside 
the region where g'{x) > 0, which is different from the par- 
allel results in the case of 7 = 4/3 (see Cao & Lou 2009). In 
Figure 5, we further exhibit a higher order g"*"— mode. It is 
clear that the oscillatory wave patterns are also mostly con- 
fined inside the collapsing stellar core, but the waves seem to 
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Figure 6. An example of g~— mode with A = 0.001, 7 = 1.30, 
I = 2 and a background profile g{x) given by expression I I39I I. 
The moving boundary of the stellar core is at xi, = 6.9561. The 
oscillations are confined inside the region where g'{x) < 0. 

overcome the barrier region (e.g. Lou 1995a,b, 1996) where 
g'{x) is negative, and are amplified when approaching the 
outer boundary of the stellar core. In the case of 7 ^ 4/3, 
however, this wave tunneling does not happen. This effect 
seems to be due to the increase of 7 value. Another feature 
of these deeply trapped g"*"— modes is that the radial oscil- 
lation dominates, while the horizontal oscillation decreases 
inside the collapsing stellar core. 

Figure 6 shows one of the g~ —modes in the case of 7 = 
1.30, as an example of the cases 7 < F = 4/3. The oscillatory 
behaviours of g~ —modes are mostly confined within a region 
of g'{x) < 0. 

In Figure 7, we show a lowest order unstable g^— mode 
for A = 0.006, 7 = 4/3, / = 1 and g{x) prescribed as 

g{x) = 1 + O.OOla;^ exp(-xV2) . (40) 

This g^— mode has an eigenvalue m = —0.00014539. Re- 
calling the eigenvalue definition of m = 3Xn(3n — l)/2, we 
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g(x)=1 +0.001 xV'^ 
y =4/3, >i=0.006 
m=-0.00014539 
1=1 



Figure 7. An unstable lowest order g"*"— mode with A = 0.006, 
7 = 4/3, Z = 1 and g{x) given by expression 1 140 I I. The contracting 
boundary of the collapsing stellar core is at xi, = 8.6598. The 
corresponding n value is either 0.3163 or 0.0170. This lowest order 
unstable g"*"— mode with 1 = 1 may describe a kick of nascent 
collapsed compact object (e.g. PNS or stellar mass black hole left 
behind in a SN remnant). 



find either n = 0.3163 or n = 0.0170. The former is larger 
than 1/6 and corresponds to an unstable temporal growth. 
As will be discussed in Section 4.3 presently, the I = 1 mode 
can lead to a displacement of the very centre of the stellar 
core. Such unstable core displacement along with nonlinear 
dynamic evolution may give rise to systematic movement of 
the collapsed core during the emergence of a reverse shock. 
Note that this mode is mostly confined inside the collapsing 
stellar core, so the displacement of the centre of the core will 
be more pronounced than that of the outer part. We there- 
fore propose that the rapid increase of this unstable central 
core displacement may be ultimately responsible for the ini- 
tial kick process of nascent PNSs or stellar mass black holes 
produced in SN remnants (see Section 4.3 and estimates in 
Cao & Lou 2009). 

Other assigned evolutions or distributions of the specific 
entropy function g{x) have also been explored, including the 
following three different forms 



O.Olx^ 
0.01a:2 I 

O.Ola;^ 



2(0.01a:2 -1^ 1) ' 
g{x) = 1 + O.lx'^ix - 2)(a; - 4) exp(-a:V9) 



(41) 

(42) 
(43) 



The results of numerical solutions for 3D general poly- 
tropic perturbations are consistent with the earlier conclu- 
sion that g^— modes are mainly confined within the regime 
with g'{x) > 0, and that g~— modes dominate in the region 
of g'{x) < 0. The specific oscillatory wave pattern may be 
infiuenced by the choice of the polytropic index 7 for 3D per- 
turbations. Regarding the last expression of g{x) as defined 
by equation (|43p or other similar expressions for which the 
first derivative g'{x) crosses the value more than once, the 
region g'{x) > [g'{x) < 0] are divided into several intervals 
in the x domain, the corresponding g+— modes (g~— modes) 
appear to have the trend to concentrate inside the inter- 
val nearest to the stellar interior core. This indicates that 



g— modes are mainly trapped deeply inside the central re- 
gion of a collapsing stellar core. 

4.3 Physical Nature of Perturbation Instabilities 

We have shown in Figures 2 and 3 that the eigenvalues of m 
for g— modes can be quite close to especially for higher ra- 
dial orders (i.e. the anti-Sturmian property). As the specific 
entropy profile g{x) tends to constant 1 and 7 approaches 
r = 4/3 simultaneously, the absolute value of buoyancy fre- 
quency squared lA/"^ | becomes smaller and smaller, squeezing 
all eigenvalues m of g— modes to converge towards zero. 

As shown in Section 4.1, we reveal a type of unstable 
perturbation solutions with eigenvalue m = (n = 1/3) 
in that very limit of g{x) = 1 and 7 = 4/3. Growing with 
time t, these are purely vortical perturbations with e^ and 
e^ components [see eq. (|25|l ]: meanwhile mass conservation 
(|37|l should be satisfied in the absence of mass density per- 
turbation (see Appendix A for details). In other words, the 
three components of vorticity perturbation (V x vi) all be- 
have in the same temporal manner of unstable growth [see 
eqns (|30p and (|31|l ] with arbitrary dependence of (x, 9, tp). 
Physically, they describe unstable vortical convections and 
circulations in a dynamically collapsing stellar core. It would 
be more convenient and precise to regard them as separate 
vorticity modes of instability. 

Complementarily, we may equally view such vortical 
perturbations as infinitely degenerate g— modes with eigen- 
value m = (n = 1/3). In reference to equation p4|l . we 
can imagine the dimensionless specific entropy g{x) to ap- 
proach constant 1 and the perturbation polytropic exponent 
7 to approach the background F = 4/3; given the fact that 
perturbation variable Wr — A{xwt) / Ax (related to ee and e^ 
components of vorticity perturbation) does not vanish for 
g— modes, it is then clear that all eigenvalues of m concen- 
trate towards 0. In that limiting regime, all g— modes would 
become completely degenerate. This also explains why vorti- 
cal perturbation solutions obtained in Section 4.1 can have 
arbitrary vortical configurations as long as that the mass 
conservation (|37p is satisfied without inducing mass den- 
sity perturbation. Therefore we may state that amplitude 
growths of g— modes in general are primarily governed by the 
angular momentum conservation (see Appendix D) . This im- 
plies that perturbation amplitudes of g— modes are roughly 
inversely proportional to the radius of the collapsing stellar 
core. This conclusion should remain valid for a stellar core 
collapse with background polytropic exponent F 7^ 4/3, for 
there should also exist amplification of perturbations result- 
ing from the angular momentum conservation. In general, 
growth rates of g— modes are also modified by the variable 
specific entropy evolution/distribution g{x) as well as the 
corresponding A value and the difference between the pertur- 
bation polytropic exponent 7 and the background polytropic 
exponent F = 4/3. Sufficiently high-order g+— modes and 
all g~— modes are unstable, that is, they grow faster than 
the temporal factor |t[~'^''® for stable acoustic p— modes and 
f— modes associated with a self-similar stellar core collapse 
(GW; Lou & Cao 2008). 

Such g— mode instabilities mainly governed by the an- 
gular momentum conservation might have a significant in- 
fiuence on the pre-SN core collapse stage and may affect 
the appearance of the PNS born during the stellar core col- 
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lapse and the emergence of rebound shock as well as an SN 
explosion. The effects of these unstable g— modes and vor- 
tical circulation modes during a stellar core collapse might 
have been omitted in previous multi-dimensional simulation 
works (e.g. Buras et al. 2006a; Dessart et al. 2006; Nordhaus 
et al. 2010b), because of their characteristic feature of weak 
or no density fluctuations. Many numerical simulations show 
that the neutrino luminosity is not high enough to trigger 
an SN explosion, ft is considered physically plausible that 
the neutrino driven convection plays a key role in an SN ex- 
plosion. Our fast growing hydrodynamic convective g-modes 
and/or vorticity modes during a stellar core collapse, if prop- 
erly included, might supply additional strengths for the PNS 
neutrino convection and circulation. These unstable pertur- 
bations will naturally destroy the spherical symmetry of the 
collapsing stellar core and give rise to subsequent aspherical 
nonlinear dynamic interactions. The resulting asymmetric 
disturbances might actually be connected to and affect the 
initiation and evolution of SASI, although SASI is primarily 
of acoustic nature and operates between PNS and the stalled 
bounce shock. In fact, g— mode oscillations of a PNS may 
also accompany SASI during the bounce phase (e.g. Burrows 
et al. 2006). Our internal g— mode, vortical and convective 
instabilities exist during the core collapse irrespective of de- 
tails of the specific entropy distribution/evolution. 

Among all these unstable g— modes, those with I — 1 
allow a non-vanishing bulk flow velocity at the very cen- 
tre of the progenitor stellar core. These unstable g— modes 
may contribute directly or indirectly to the observed high 
peculiar speed of a new-born PNS (remnant of a stellar 
core collapse), which may reach values of several hundred 
to more than one thousand kilometers per second. As an 
estimation, we assume the initial gas bulk flow velocity of 
the stellar core relative to the center of mass of the entire 
star to be ~ 200 m , and the radius of the collapsing 
PNS contracts by a factor of ~ 10~^, then the stellar core 
will gain a velocity around ~ 200 km s~^ as a result of the 
rapid growth of I — 1 g— mode perturbations, provided that 
the linear approximation still remains valid|j The smaller 
this contraction factor, the higher the PNS kick speed; the 
phase of PNS 1 = 1 g— mode oscillation, the proper timing 
of stellar core collapse and the overall flow environment to- 
gether should be also important to determine the net kick 
efficiency. The initiation of such a collapsing PNS kick may 
actually occur prior to the emergence of a bounce shock. 
Recent numerical simulation studies suggest that the neu- 
tron star kick may be produced by the gravitational pull of 
the post-bounce anisotropic slow-moving ejecta (see figure 7 
of Scheck et al. 2006) which may result from the neutrino- 
driven convection (e.g. Scheck et al. 2004; Nordhaus et al. 
2010a) i For example, based on extensive 2D model simu- 
lations with approximate neutrino transport and boundary 
conditions that parameterize the effects of the contracting 
NS and strong enough neutrino heating for neutrino-driven 

^ For a collapsed stellar core temperature of several times 10^^ K 
(e.g. Bethe 1990), the PNS core sound speed would be in the 
order of ~ 3 X 10* km s~^ and the estimate based on the linear 
approximation here should be justifiable. 

^ Parts of the work by Scheck et al. (2006) have already been 
presented in Scheck et al. (2004), but a detailed description of 
both their methods and results are given in the former. 



SN explosions, Scheck et al. (2006) compared their results 
with those of Burrows et al. (2006). Started 20 ms after the 
core bounce, the simulation results of Scheck et al. (2006) 
can lead to a dominance of dipole {I = 1) and quadruple 
(/ = 2) modes in the explosion ejecta, provided that the on- 
set of the SN explosion is sufficiently slower than the growth 
time scale of the low-mode instability; they show that global 
anisotropies and large NS kicks can be obtained naturally 
in the framework of the neutrino-driven SN explosion mech- 
anism due to the symmetry breaking by non-radial hydro- 
dynamic instabilities, without resorting to rapid rotation, 
large pre-coUapse perturbations in the iron core, strong mag- 
netic fields, anisotropic neutrino emission associated with 
exotic neutrino properties, or jets etc. As a complement to 
the work of Scheck et al. (2006), Nordhaus et al. (2010a) 
performed a 2D axisymmetric radiation-hydrodynamic sim- 
ulation of the collapse of a non-rotating 15Mq progenitor 
core and can naturally capture the PNS formation and sub- 
sequent off-center acceleration during a delayed, neutrino- 
driven, anisotropic SN explosion (their SN explosion is ar- 
tificially induced by including additional neutrino luminos- 
ity to the calculation); at the end of their simulation, the 
PNS has achieved a velocity of ~ 150 km s~^ and is still 
accelerating at ~ 350 km s~^. In the simulation of Burrows 
et al. (2006), the acoustic energy input from PNS g-mode 
oscillations was crucial for the explosion of an IIMq pro- 
genitor. Nonspherical accretion was found to give rise to 
core g-mode oscillations at late times ( > 300 — 500ms) af- 
ter core bounce, providing a significant amount of acoustic 
power for SN shock. Actually, g-mode oscillations are also 
present in the outer layer of a NS according to Scheck et 
al. (2006), but the amplitudes are modest and do not lead 
to the strong effects of Burrows et al. (2006). It is possible 
that Scheck et al. (2006) simulations might underestimate g- 
mode effects which would require the inclusion of the entire 
neutron star without excising the central core and thus the 
dynamic coupling among accretion, core motion, and deep g- 
mode generation in a self-consistent manner. According to 
our estimates above, the core collapse hydrodynamic con- 
vective instabilities may also have a significant impact on 
the development of the anisotropic ejecta and the asymmet- 
ric bounce shock, and may also contribute directly to the 
PNS kick, providing a pre-existing background velocity dur- 
ing the formation of a nascent PNS. It would be desirable 
that these sources of g— mode instabilities during core col- 
lapse be further explored by 3D numerical simulations. In 
addition, those unstable g— modes with I = 2 might induce 
a collapsing stellar core to break up into two pieces, form- 
ing binary PNSs or other combinations around the centre 
of the pre-SN progenitor, while g— modes with higher I and 
under favorable conditions may violently destroy the stellar 
core, shredding it into many smaller pieces and thus lead- 
ing to the destruction of a PNSllI The possibility of such a 
situation might be low, as several multi-dimensional numer- 
ical simulations show that low— ? modes appear most likely 
dominant. 



^ Although highly speculative, the absence of evidence over more 
than two decades for a central compact object in the remnant of 
SN1987A might be indicative of such a possibility of broken core. 
Admittedly, there could be other alternatives to resolve this issue. 
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Based on our model analysis, the vorticity modes of in- 
stability appear unavoidable during a stellar core collapse; 
this would happen irrespective of the specific entropy evolu- 
tion/distribution g{x) including the conventional polytropic 
case of g{x) = 1. Depending on sources of seed fluctuations, 
rapid core spin, circulations of various scales, and vortical 
convections would go along with the stellar core collapse 
in a generic manner. Physically, such 3D vortical motions 
can naturally lead to core convective turbulence and thus 
MHD dynamo actions for producing intense magnetic fields 
in PNSs and in stellar mass black holes. Conceptually, we 
may specifically separate the stellar core rotation from the 
rest of radial vorticity perturbation [see eq. (|3ip ] and envi- 
sion the scenario of nonlinear Rossby waves (e.g. Haurwitz 
1940; Lou 1987, 2001). 



5 SUMMARY AND CONCLUSIONS 

3D general polytropic perturbations in a F = 4/3 self- 
similar collapsing stellar core have been examined. We com- 
bine two recent generalized formulations adopted by Cao 

6 Lou (2009, 2010) together to explore even more gen- 
eral cases of 3D perturbations. The specific entropy evolu- 
tion/distribution profile g{x) is allowed to vary from the stel- 
lar centre to its outer contracting boundary. The special sub- 
case of uniform g{x) = 1 corresponds to a conventional poly- 
tropic gas sphere (GW; Lou & Cao 2008). The polytropic 
exponent 7 of 3D perturbations is generally allowed to be 
different from the background polytropic exponent F = 4/3. 
Not surprisingly, these general conditions give a non-zero 
Brunt- Vaisala buoyancy frequency squared Af^ as defined by 
equation (|38|l , and allow for two groups of g— mode perturba- 
tions to occur, viz. g"""— modes and g~— modes under proper 
situations. Meanwhile, acoustic p— modes and f— modes can 
also be clearly identified from the results of numerical cal- 
culations; again, they remain stable. 

We first revealed a class of unstable radial vorticity per- 
turbation solutions that are incompressible and can be sep- 
arated out independently from the rest of compressible per- 
turbation variables under fairly general situations, i.e. irre- 
spective of the specific entropy evolution/distribution profile 
g{x). Such unstable vortical perturbations include the spin 
of the collapsed stellar core and all kinds of sufficiently slow 
vortex motions tangent to spherical layers in the stellar core. 
Their velocity amplitudes are inversely proportional to the 
radius of the collapsed stellar core, as a result of the angu- 
lar momentum conservation. This appears to be a natural 
mechanism for producing fast spins of PNSs or stellar mass 
black holes in the eventual nonlinear evolution. Conceptu- 
ally, one may further separate the core spin and other radial 
vorticity perturbations; their mutual nonlinear interactions 
and evolution may be perceived. We therefore suggest that 
Rossby waves may occur and be amplified in nonlinear dy- 
namic processes. Such Rossby waves may be also connected 
with the spin down process of a nascent neutron star. If the 
collapsing stellar core is magnetized, significant MHD dy- 
namo actions may also take place to amplify magnetic fields 
during the core collapse process involving considerable tur- 
bulent convections and circulations. 

By the inverse iteration method (e.g. Wilkinson 1965), 
we then solve the remaining general polytropic compressible 



perturbation equations (|32|) — (I34|) numerically with speci- 
fied boundary conditions (I35|l . and examine the eigenval- 
ues and eigenfunctions of perturbations, including acous- 
tic p— modes, surface f— modes and internal g— modes, 
which can be further divided into two subclasses: internal 
g"*"— modes and internal g~— modes. Sufficiently high radial 
order g'''— modes and all g~— modes show instabilities dur- 
ing the dynamic core collapse process, leading to convective 
motions of circulation within the collapsing stellar core. In 
contrast, acoustic p— modes, f— modes and sufficiently low 
order g^— modes remain stable. For a certain specific en- 
tropy evolution/distribution, there may exist g— modes that 
are drastically oscillatory deeply trapped in the collapsing 
stellar core. 

The mechanism of unstable growths of g— modes is iden- 
tified. We realize that the fast temporal increase in am- 
plitudes of g— modes are mainly governed by the angular 
momentum conservation. It is easy to understand this in- 
tuitively, for g— mode perturbations are mainly convective 
motions inside the dynamically collapsing stellar core, in 
which the gas runs in circulations. During the core collapse, 
the velocity of the gas fiow will increase accordingly (see 
Appendix D for more details) . In a more realistic model for 
stellar core collapse, both k and F may vary in time. How- 
ever, during a long period before the end of the collapse of 
stellar core, numerical simulations (e.g. Burrows & Lattimer 
1986; Van Riper & Lattimer 1981; Bethe 1990) show that 
the variations of k and F are slow and not obvious, and 
our conclusions may still be valid approximately. The am- 
plitude of g— mode perturbations are more or less inversely 
proportional to the radius of collapsing stellar core. When 
the specific entropy distributes uniformly and the pertur- 
bation polytropic exponent 7 = F = 4/3, this temporal 
increasing factor of perturbation amplitude becomes exact, 
and all the g— modes become completely degenerate in this 
regard. Alternatively, we can regard these unstable modes as 
3D vorticity perturbations constrained by the mass conser- 
vation. They may bear physical consequences for the central 
compact object, be it PNSs or stellar mass black holes. 

Given the idealizations and approximations of our 
model, these unstable g— modes may offer certain valuable 
clues for SN model simulations. They affect the geometry of 
the rebound shock wave emerged after the core bounce be- 
cause of asymmetries and enhance the neutrino driven PNS 
convection. Most simulation works do not explicitly include 
these convective g-mode perturbations (e.g. Dessart et al. 
2006; Nordhaus et al. 2010b). Buras et al. (2006b) added in 
their simulation an arbitrary density perturbation before the 
core collapse, and do not observe significant difference. The 
density perturbations, however, are likely more acoustic and 
grow slowly, because the convective g-modes involve weak 
density fiuctuations. Among the various g-modes, the 1 = 1 
g— modes with central movement may possibly be respon- 
sible for the origin of the high kick speed that a new-born 
neutron star acquires, or serve as the source of producing 
anisotropic stellar ejecta which may possibly gravitationally 
accelerate a nascent PNS (Scheck et al. 2006; Nordhaus et 
al. 2010a). Nonlinear evolution of those unstable g— modes 
with I higher than 1 may possibly break a collapsing core 
into two or multiple pieces, forming a binary system of neu- 
tron stars or more fragments as to even ultimately prevent 
the formation of a neutron star. 
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Instabilities of radial component vorticity perturbation 
can lead to rapid rotations of central compact objects, 
Rossby waves, and circulations over spherical surfaces in- 
side a core. Instabilities of 3D vorticity perturbation when 
g{x) = 1 and 7 = F = 4/3 can give rise to convective tur- 
bulence in the core and MHD dynamo actions to sustain 
violent magnetic activities. 
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APPENDIX A: EQUATION (35) IS NECESSARY 
GIVEN 7 = 4/3, G{X) = 1 AND M = 

In the special case of 7 = 4/3, g{x) — 1 and m — for a 
conventional poly tropic gas, equation (|34l) is automatically 
satisfied. In the first subcase of n = 0, equation (|37|) comes 
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out automatically from equation (|32|) . When we focus on the 
other subcase of n = 1/3 satisfying the same requirement 
m = 0, we derive from equation (|33p the following relation 



/ d 2 
ax 



l{l + l)f 



-Wt 



Sf'wr 



meanwhile equation (|32|) yields 



dx^ 
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{x^pi 
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In the interval {xt, +oo), we define f{x) = 0. To avoid sin- 
gularity, we naturally require xipi to remain finite at a; = 
and X +00. Equation (|32} should be regarded as a Sturm- 
Liouville problem within the interval [0, +oo) for determin- 
ing the eigenvalues l{l + l). The eigenvalues / are determined 
by the specific profile of the dynamic background density 
function f{x) and in general cannot be exactly equal to an 
integer. This means that in general we could only have the 
trivial solution xipi = for integral values of I. Another ar- 
gument proving the necessity of xtpi = can be made in 
the following for those cases I ^ 2. Since xipi should remain 
finite at a: = to avoid singularity of a proper perturba- 
tion solution, xipi oc x should hold as x ^ 0+. We also 
need to require xipi (x x~' when a:: — >■ -|-cxd so that ipi ap- 
proaches zero at infinity. This means that the algebraic sign 
of d{xipi)/dx should necessarily reverse when x varies from 
0"*" to -|-oo. Suppose that xipi is positive when a; — > 0"*", then 
its first derivative is also positive. Ifowever, from the nu- 
merical results of f{x) it is found that for I > 2, inequality 
l{l + l)/x^ — 3/^ > holds in the entire interval (0, -f-oo), 
which indicates that d^ {xipi)/dx'^ remains positive as long 
as xipi is positive. There seems no chance for the sign of 
d(a;V'i)/da; to reverse. Then xipi must vanish. In conclusion, 
we show that xipi = 0, and according to equation (Al), the 
necessity of equation p7p is warranted. 



APPENDIX B: ORTHOGONAL SET OF 
EIGENFUNCTIONS WITH DISTINCT 
EIGENVALUES M 

Perturbation eigenfunctions with different eigenvalues of 
m = 3An(3n — l)/2 can be proven to be mutually orthogo- 
nal to each other. In the proof below, we use superscripts 
and to distinguish different eigenfunctions and eigenval- 
ues. From perturbation equations (|22p and ([23]), we readily 
find 

V'5^=2V (/Vi) ; (Bl) 
or 

a substitution of equations (|23|) and (|24|l into equation (|21[) 
then yields 



4 .3 
3m/ V. 



7»/'(V-vi)+vi.V(f;/* 



;V.v.+v..^)v(5/) + 2/V^ 



(B2) 



We now evaluate the spatial integral below with the help of 
equations (|B1|I and (IB2|l . namely 



4 (1) / r3 (1) 



Vl 



(2) 



dV 



. (1) 



75/ V • Vl' * + Vl 
V-vi^') -hvi 

dipi'^ 

dr 



(1) 



V(<7,/1 



V-vi'^Mv" 



3V/ 



V(ffr).vi'^MF 



-V- [Ai'^'jdl/ 



. (2)^ 



j |7ff/[V.viW](V.vi( 

+ [vi(^'-V(3.f*)](V.vi(^)) 
+ [vi(^'-V(3/)](V.viW) 

(2)1 3V/ ■ V(ff/") 



+ [er ■ vi'^'j [er ■ Vl 



(B3) 



where dV is the volume element for the spatial integration. 
Noting that the RHS of equation (|B3|l is symmetric with 
respect to two superscripts '^^^ and we have 



[^a)_^{2)] j /Vi'^'-vi(^'di/ = o 



(B4) 



and the orthogonality of eigenfunctions with different eigen- 
values of m is therefore proven. 



APPENDIX C: INVERSE ITERATION SCHEME 
FOR SOLVING LINEAR ODES (30)-(32) 

We give here a procedure description of the inverse iteration 
scheme (e.g. Wilkinson 1965) as adapted to our numerical 
computations of 3D general polytropic perturbations in a 
self-similar dynamic stellar core collapse of spherical sym- 
metry and F = 4/3. 

First, we specify a chosen form of the specific entropy 
function ^(a::) and employ the fourth-order Runge— Kutta 
scheme (e.g. Press et al. 1996) to integrate background ODE 
P5|l for f{x) from x = and thus obtain the dynamic 
background function f{x) with the initial or boundary con- 
ditions /(O) = 1 and /'(O) = 0; this f{x) is closely related 
to the background mass density profile. For F = 4/3, the 
contracting boundary of the stellar core xi, is determined by 
the realizable requirement f{xb) = in a parameter range of 
< A ^ Ac for a proper A value. Our goal is to reliably solve 
linear DDEs (|32|l ~ (|34p for 3D general polytropic compress- 
ible perturbations within the interval [0, a:;,] by determining 
the eigenvalues of m = 3An(3n— l)/2 and the corresponding 
perturbation eigenfunctions for given g{x), f{x), A, xi,, I and 
7 7^ F = 4/3 in general. 

In order to implement the inverse iteration scheme, we 
need to discretize linear DDEs (|32p — p4p with a proper mesh 
in the interval (0, a;^). More specifically, we first divide the 
entire interval evenly into A'' small intervals and define the 
interval size Aa: = h = Xt/N; we can then construct a vec- 
tor w of length 2{N ~ 1) for Wr and wt arranged together 
contiguously in order, namely 



Wr{h),Wr{2h), ■ ■ ■ Wr{{N - l)h), 
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Wt{h),Wt{2h),---wt{{N ^l)h) 



(CI) 



We also need to cast every differential operator into a 
matrix form. The matrix for a function serving as coefficients 
in the ODEs, say y{x), would then take the form of 



Y = Diag y{h),y{2h),- ■ ■ y{{N ~ l)h) 



(C2) 



which is a matrix of dimension (A^ — 1). By notation ^Diag\ 
we refer to a diagonal matrix. The matrices D and for 
the differential operators d/dx and d? /dx^ are respectively 



D 



and 



\ 



-1 ij 



/ -2 1 \ 
1 -2 1 

1-2 1 

V 1 -2 / 



(C4) 



both are of dimension (A'^ — 1). For different boundary condi- 
tions, these matrices above may need modifications accord- 
ingly. Note that ODE ^ can be substituted into ODE 
to eliminate the dependent variable nipi. With these pro- 
cedures, we could finally achieve the discretization scheme 
of linear ODEs (P^ — (Pl|l . The resulting finite diflFerence 
scheme can be simply cast into the following form of 



Ai A2 

As Ai 



A-w 



(C5) 



in which Ai {i = 1, 2, 3, 4) are four matrices all of dimension 
(A*' — 1) that do not involve the eigenvalue parameter m; 
here, eigenvalue parameter m appears only on the RHS and 
is to be determined for discrete eigenvalues. 

We can then apply the inverse iteration method to solve 
this eigenvalue problem. We begin to specify an arbitrary 
vector Wo and guess an approximate eigenvalue mo, and 
then rewrite the matrix equation (IC5|) as 



{A — mo)w = (m — mo)w 



(C6) 



Suppose that all the eigenvalues and eigenfunctions of ma- 
trix equation (ICSP are respectively rrii and Wi (i = 1, 2, ■ ■ ■), 
among which the value of rrij is closest to the initially cho- 
sen value mo. Then, wo can be expanded as wo = ^ . OiW;, 
where Ui stands for a set of coefficients, and it is thus easy 
to demonstrate 



APPENDIX D: MAGNITUDE GROWTH OF 
G-MODES DURING A CORE COLLAPSE BY 
ANGULAR MOMENTUM CONSERVATION 

The angular momentum in a certain spatial volume V of 
first-order precision can be cast in the integral form of 

= y p(r X u)d^r = aa J /^(x x vi)d^x 

= Ae''^'' 1 /^(x X vi)d^x , (Dl) 
Jv 

where ^ is a constant. Taking the derivative of the angular 
momentum in terms of the logarithmic time r = — In |t| and 
using equation (|2ip . we obtain 



(C3) ^ 
dr 



Ae 



/3 



t/3 



d 1 



d^x 



(D2) 



From equations (|22|l — p4|l . it immediately follows that 



V'^Vi = 2 / V • (/'vi)dr 



/3i 



7V-(/^vi)-hvi- (^(4-37) 



f 9 ) 



(D3) 



dr.(D4) 



Since V'i/'i is a potential vector field, integral (D3) can be 
further simplified as 



ViAi = 2 / if\i)pdr 



(D5) 



where (/'^vi)p stands for the potential flow part of the vector 
fleld /^vi. Since for g— modes the vector field /^vi (which 
is proportional to the momentum density) is close to diver- 
gence free [which becomes exact when 17(2;) = 1 and 7 = 4/3, 
see equation (|37p ]. Vi/'i of g-modes should be quite small. 
For nearly constant g{x) and 7 ~ 4/3, it is clear that /3i 
is also a small quantity. We therefore have an approximate 
angular momentum conservation law 

dJv 







(D6) 



which indicates (d/dr — l/3)vi « 0. It is thus the conserva- 
tion of angular momentum that mainly governs the temporal 
increase of the velocity perturbation amplitude. 

This paper has been typeset from a Tff)i/ ETJ5X file prepared 
by the author. 



{A — mo) "wq — >■ aj{mj ~ mo) "wj as n — > 00 . (C7) 

Therefore, by iterating an enough number of times, we can 
select the eigenfunction Wj and determine the corresponding 
eigenvalue rrij = {rrij — mo) -I- mo. Changing the value of 
mo, we would be able to select other eigenfunctions and 
eigenvalues of the perturbation problem. 



